Influence of the Ground-State Topology on the Domain- Wall Energy in the Edwards-Anderson ±7 

Spin Glass Model 



O 
O 

(N 

O 

Q 
c 



F. Roma and S. Risau-Gusman 

Centro Atomico Bariloche and CONICET, R8402AGP San Carlos de Bariloche, Rio Negro, Argentina 

A.J. Ramirez-Pastor and F. Nieto 

Departamento de Fi'sica, Universidad Nacional de San Luis and CONICET, Chacabuco 917, D5700BWS San Luis, Argentina 

E.E. Vogel 

Departamento de Fi'sica, Universidad de La Frontera, Casilla 54-D, Temuco, Chile 
(Dated: February 6, 2008) 

We study the phase stability of the Edwards- Anderson spin-glass model by analyzing the domain-wall energy. 
For the bimodal ±7 distribution of bonds, a topological analysis of the ground state allows us to separate the 
system into two regions: the backbone and its environment. We find that the distributions of domain-wall 
energies are very different in these two regions for the three dimensional (3D) case. Although the backbone 
turns out to have a very high phase stability, the combined effect of these excitations and correlations produces 
the low global stability displayed by the system as a whole. On the other hand, in two dimensions (2D) we 
find that the surface of the excitations avoids the backbone. Our results confirm that a narrow connection exists 
between the phase stability of the system and the internal structure of the ground-state. In addition, for both 3D 
and 2D we are able to obtain the fractal dimension of the domain wall by direct means. 
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The spin glass state has been studied extensively during the 
last thirty years, but the role played by low energy excitations 
is still a matter of debate. These excitations are crucial to un- 
derstand the nature of the ordering of the spin glass phase. 
Most studies have focused on the predictions of two theo- 
ries: the replica-symmetry breaking (RSB) picture [1] and the 
droplet picture |2y. RSB, rigorously true for the Sherrington- 
Kirkpatrick model of spin glasses, predicts that there are ex- 
citations which involve flipping a finite fraction of the spins 
and, in the thermodynamic limit, cost only a finite amount of 
energy. The fractal dimension of the surface of these excita- 
tions, d s , is expected to be equal to the space dimension, d. 
On the other hand, in the droplet picture the lowest energy ex- 
citations of length L have d s < d and typically cost an energy 
of order L e (9 is known as stiffness exponent). Thus, contrary 
to RSB, the droplet picture predicts that excitations involving 
a finite fraction of spins cost an infinite amount of energy in 
the thermodynamic limit. 

The exponent 9 plays a central role in this debate. It is usu- 
ally calculated by using the concept of defect energy, AE = 
E" — E 1 ', which is the difference between the ground-state 
(GS) energies for antiperiodic (E") and periodic (E p ) bound- 
ary conditions, in one of the directions of a af-dimensional 
system of linear size L. In ferromagnetic systems, AE ~ L e , 
with 9 = d s = d — 1, because the induced defect is a (d — 1)- 
dimensional domain wall with all its bonds frustrated. For 
spin glasses, the average over the distribution of bonds (de- 
noted by [...]) must be taken and the scaling ansatz becomes 
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Assuming that, because of frustration, the defect energy is the 
sum of many correlated terms of different signs, Fisher and 
Huse |01 have shown that for spin glasses 9 < 

It is well known that the Edwards-Anderson (EA) ±7 



modelyl] has a degenerate GS. In addition, it has been shown 
|H|5l0that, for each realization of the disorder in two dimen- 
sions (2D), there are bonds which are either always satisfied or 
always frustrated in all the GSs. These bonds define the Rigid 
Lattice (RL), or backbone of the system. Spins connected by 
it are called solidary spins (the rest are denoted non-solidary). 
Numerical studies in three dimensions (3D) of a similar struc- 
ture, the Diluted Lattice [J7[| (defined by the satisfied bonds of 
the RL), seem to confirm the existence of the RL in the ther- 
modynamic limit. 

Recently, an interesting connection has been found be- 
tween the RL and both thermodynamic and dynamic prop- 
erties fjLls HlOll . For example, in Ref. ifioll the slow and fast 
degrees of freedom associated to the out-of-equilibrium dy- 
namics of the 2D EA ±7 model were shown to be closely 
related to the solidary and non-solidary spins, respectively. In 
this context, one expects that the contributions of the back- 
bone and its environment on the domain-wall energy will be 
rather different. The present letter represents a step in that 
direction. 

In this work we show that in 3D a significant portion of 
the domain wall is contained in the backbone. Moreover, in 
this region we find that the defect energy is positive and has 
a strong size dependence, similar to what happens in a ferro- 
magnetic system. The rest of the system has a negative de- 
fect energy and could therefore be considered as an excited 
phase. The sum of these two defect energies results in the 
cancellation effect responsible for the weak size dependence 
observed for the domain-wall energy fill [l2il . On the other 
hand, we find that in 2D the portion of the domain wall inside 
the backbone is vanishing. We have also obtained the fractal 
dimensions of domain wall for both 2D and 3D. 

We start by considering the Hamiltonian of the EA model 
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TABLE I: Parameters of the simulation for each lattice size L in 3D. 
Nj is the number of temperatures used in the parallel tempering, cho- 
sen inside the interval T m - m to T max (temperatures in units of J /kg), 
MCS is the number of Monte Carlo steps needed to reach the GS and 
Ns is the number of samples used. Ngs is the mean number of GSs 
used to approximate the average ((...)gs) f° r the antiperiodic system. 
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for spin glasses on square and cubic lattices, 

ff=J^/ya / a /l (2) 

(>',;) 

where a,- = ±1 is the spin variable and indicates a sum 
over nearest neighbors. The coupling constants are indepen- 
dent random variables chosen from a ±J bimodal distribution. 

For 3D lattices and periodic boundary conditions, we de- 
termine the RL by using an improvement of the algorithm 
introduced in Ref. 11311 . where parallel tempering 111411 has 
been implemented for reaching the GS. The present scheme, 
called Rigid Lattice Searching Algorithm (RLSA), allows to 
obtain true GSs up to L = 12, where L is the lattice size. How- 
ever, the maximum size used in this work is L — 8, because 
to obtain the RL, the GS must be reached 3N (N = L d ) times. 
An important point is that the set of parameters used in the 
parallel tempering (see Table I) has not been chosen to equili- 
brate the system, but to reach quickly a GS configuration. To 
check this, for each lattice size we have compared the average 
GS energy per spin, eo, calculated from our algorithm, with 
the value reported in the literature lfl5l [l6ll . For example, for 
L = 8 we obtain eo = —1-780(1), in good agreement with the 
value e = -1.7802(5) reported in Ref. fm . 

For 2D lattices, we have used a different algorithm. It is 
well known |[T5ll that the problem of finding the GS for a 
2D lattice with at least one free boundary condition can be 
mapped to a minimum weighted perfect matching problem, 
for which very efficient algorithms exist. To implement the 
RLSA we have used one of these routines, which has allowed 
us to calculate the RL up to L = 100. 

As mentioned above, the RL is formed by bonds (rigid 
bonds) which are either always satisfied or always frustrated 
in all the GSs. The remaining bonds, called flexible bonds, 
form the flexible lattice (FL). This allows us to write the 
Hamiltonian ©as: H = H r +Hf. Subindex r (/) refers to 
the Hamiltonian restricted to only rigid (flexible) bonds. 

To calculate the defect energy we write the GS energy of a 
particular sample (p), as 

E P =E P +E P . (3) 

Note that, although the GS is degenerated, E p , E? and E p 
remain constants on all configurations of the GS. Next, a new 
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FIG. 1: Distribution P, (full line) and P f (dotted line) for 3D. The 
energy is given in units of J. The inset shows the scaling of [AE r ], 
[AE r ] — c and [AE m ] — c. Error bars are smaller than the symbols. 

sample is generated (a), by introducing antiperiodic boundary 
conditions in one direction. Its energy can be written as 

E a = {E?) GS +(Ef) GS , (4) 

where the subindexes r and / correspond to the restriction of 
the Hamiltonian to the bonds that form the RL and FL of the 
periodic system, and (...)gs denotes an average over all the 
GSs of the antiperiodic system 11711 . Notice that this average 
is necessary, because now E" and E" are not constants in all 
the GSs. Using Eqs. (01 and we write the defect energy as 

AE = AE r + AE f , (5) 

where AE r = (Ef) GS - E p and AE f = {Ef) GS - E p f . 

In 3D, we have measured [lAZsl] for each lattice size (due to 
symmetry arguments, [AE] = 0) and we have obtained a finite- 
size dependence with a small stiffness exponent 9 w 0.2, in 
agreement with the result of Refs. ITll [l2tl . This is seen as 
an evidence of the existence of a finite critical temperature. 
We have also measured the distribution of |A£|. For the sizes 
considered, this function extends up to |AE| = 12/ (the val- 
ues of the defect energy are multiple of 47). In contrast, for 
the 3D ferromagnetic Ising model, the defect energy takes the 
value 1287 for L = 8. The small values obtained for the EA 
model indicate the presence of large sets of frustrated or sat- 
isfied bonds that compensate each other. 

The main point of this article is that these sets can be very 
clearly related to the topology of the system. This can be seen 
on Fig. Q] which shows the distributions of the contributions 
of the RL, P r (AE r ), and the FL, Pf(AEf), to the domain- wall 
energy. Two features stand out. On the one hand, AE r is al- 
most always positive. Only a small fraction of samples have 
AE, < 0, and this fraction decreases with L. Our measure- 
ments indicate that on the portion of the domain wall that 
crosses the RL, the fraction of excitations (rigid bonds always 
satisfied that appear frustrated in the GS of the antiperiodic 
system) seems to tend to 0.72. On the other hand, P r is a broad 
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distribution and extends up to high values of AE r . In addition, 
our simulations show that the distribution of AE r seems to be 
very broad also for fixed values of |AE|. For example, for 
L = 8, we have found samples with AE r = —AEf = 123.597 
or AE r = —AEf =11 .637, both corresponding to AE = 0. 

The inset in Fig. Q] shows the scaling of [AE,] with L. It 
tends towards a power law behavior of the form 

[AE,.]~L e <\ (6) 

To determine 9 r , we have fitted the data with the function 
[AE,] = c + bL 6r (which is the simplest correction to scaling). 
We obtain a good fit for c = 3/2 and 6, = 2.59(2). The inset 
in Fig. Q] shows the scaling of [AE,] — c with L. Notice that 
the exponent 9 r is an order of magnitude bigger than the usual 
stiffness exponent 9 (« 0.2 for 3D). 

The behavior of PAAEA is very similar to the one observed 
in P r , but now AEf is always negative (see Fig.[T|i. The vanish- 
ing of the average domain- wall energy implies that — [AEf] 
follows a power law with the same exponent as [AE,]. 

The previous results indicate that the proposal of separat- 
ing the system in two regions is not trivial: the sections of the 
domain-wall energy with positive (negative) sign, prevail on 
RL (FL). However, the system has a small but positive stiff- 
ness 6, because on average the defect energy [AE,] overcomes 
the defect energy [AEf] . 

Our numerical results also allow us to infer that the expo- 
nent r is equal to the fractal dimension of domain wall, d s . 
To justify this conjecture, let us consider the topological char- 
acteristics of the RL. Our simulations up to L = 8 indicate that 
the RL consists mainly of a compact percolation cluster (PC) 
Hill . The fraction of the RL that corresponds to this cluster 
converges to 0.78. The fact that the distributions P r are vanish- 
ing for negative AE r shows that compensation effects are not 
important in the RL. Assuming that the same happens in the 
PC and that it behaves as a natural box containing (on average) 
always the same fraction of domain wall as the whole RL, it 
is natural to conjecture that the area of the domain wall inside 
the RL follows a power law L* . Therefore, d s « 9 r = 2.59(2), 
which agrees with the values reported in the literature for the 
EA model with Gaussian distributed couplings lfl9ll . 

To refine this, measurements were carried out separating 
the system into two new regions: the PC and its environment 
(now, FL plus small clusters of RL). Then, the defect energy 
can be written as AE = AE pc + AE e , where subindexes pc and 
e refer to the defect energy of PC and its environment, respec- 
tively. If we assume again a power law, [AE pe ] ~ L e i' c , we 
obtain pc = 2.57(2), which gives d s = 9 pc = 2.57(2), consis- 
tent with the value quoted above. The inset in Fig. [1] shows 
the scaling of ( [AE /X ] — c) with L for c = 3/2. 

The picture in 2D is very different. The main problem is 
that, even though the RL spans a significant portion of the 
square lattice, it does not percolate l6l ll8ll . It consists instead 
of a large number of relatively small islands. The analysis of 
the difference between periodic and antiperiodic systems re- 
veals that the number of bonds in the RL that belong to the 
domain wall grows as ~ L° 6 . But the size of the domain wall 
is necessarily larger than L, which implies that the fraction of 
it that crosses the RL is vanishing with L. Geometrically, what 




FIG. 2: Matchings corresponding to the ground states of a 2D sample 
with periodic (left) and antiperiodic (right) boundary conditions on 
the horizontal direction. Points represent frustrated plaquettes. The 
grey line crossing the sample is the resulting domain wall. 



is happening is that the domain wall crosses the sample avoid- 
ing the RL islands (notice that this is further evidence that the 
RL does not percolate). Consequently, in the following we 
shall use a different strategy for studying the influence of the 
backbone on the domain-wall energy in the 2D case. 

Thus, for each sample, the domain wall depends on the pair 
of GSs (one GS for the periodic and another for the antiperi- 
odic system) that are being compared. However, with the hope 
of capturing its relevant properties, for each sample we have 
picked a single random pair of GSs. But even for a pair of 
fixed GSs the determination of the domain wall is not trivial. 
It can be shown [18] that there are some sets of bonds of the 
periodic system that appear flipped (in the sense that their sat- 
isfied bonds become frustrated and viceversa) in the antiperi- 
odic system but do not contribute to the energy change. Notice 
that this implies that they can also appear flipped in other GSs 
of the periodic system. Therefore, one must weed out these 
sets of bonds to get the correct domain wall. The plaquette 
picture [20] provides the best framework for this. 

The 2D lattice can be seen as a set of squares, called plaque- 
ttes, bounded by four bonds. If an odd number of these bonds 
are frustrated, the plaquette is called frustrated. It has been 
shown that plaquettes are frustrated if and only if there is an 
odd number of ferromagnetic bonds in the boundary I20I I2TI1 . 
Every spin configuration can be mapped to a perfect match- 
ing, which is a set of paths, made of frustrated bonds, that 
join pairs of frustrated plaquettes. The energy of the system 
is proportional to the total length (also called weight) of these 
paths. Thus, finding a GS is equivalent to finding a minimum 
weighted perfect matching. 

To define the domain wall we use, as an example, the sys- 
tem shown in Figure [2] (generalization is straightforward). To 
account for the free boundaries two plaquettes must be added, 
at the top and the bottom (outside the lattice in Fig. 13 11511 . 
The top (bottom) plaquette is called frustrated if there is an 
odd number of frustrated bonds in the top (bottom) boundary. 
The antiperiodic boundary condition changes only the frus- 
tration of these two plaquettes, which in turn leads to a new 
matching as GS. From the comparison of the two GSs a set 
of contours can be defined, where each contour is formed by 
alternating paths from each matching (see Fig. |2). All these 
contours will be closed (as the one joining plaquettes 8 to 11 
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FIG. 3: Scaling of ratio [/,-//] with 1/L for 2D. The inset shows the 
scaling of [/] with L. Error bars are smaller than the symbols. Av- 
erages have been taken over 1000 samples of each lattice size up to 
L = 100. 

in Fig. [2), except for one that runs from the top plaquette to 
the bottom one. This last contour is the domain wall (it can be 
proved that the loops do not contribute to the energy change). 

The results of simulations performed for several sample 
sizes are shown in Fig. [3] where both the average length of the 
wall, [I], and the average fraction of RL bonds in the domain 
wall, [l r /l], are shown. As anticipated, this fraction vanishes 
for large L although a very small, but nonvanishing value, can- 
not be ruled out. The average length of the domain wall fol- 
lows a scaling ([/] ~ L ds ) that allows us to find its fractal di- 
mension, with d s = 1.30(1). Remarkably, this value coincides 
with the fractal dimension reported [ 22] for the 2D EA model 
with a Gaussian distribution of bonds. 

In summary, we have studied the relevance of characteriz- 
ing the domain wall of the EA ±J spin glass by using topo- 
logical information of the GS. The defect energy on the RL in 
3D shows a behavior typical of a highly stable phase (similar 



to the 3D ferromagnetic Ising model, but with d s > 2). On 
the other hand, the FL shows a very different behavior, like a 
system in an excited state. The total defect energy AE is the 
result of these competitive and correlated effects (the whole 
system shows a low stability with a small stiffness exponent). 

The 2D case is very different. We have shown that the 
defect energy avoids the RL, lying almost completely on FL 
(which percolates in 2D). If we assume that the RL is the only 
structure able to support a stable phase, then this behavior is 
compatible with an unstable phase and a zero critical temper- 
ature. This agrees with most recent studies (see e.g. ll23ll ). 

In addition, for both dimensionalities we obtain the fractal 
dimension of the domain wall by direct means. These values 
are in good agreement with the ones reported in the literature 
for the EA model with Gaussian distributed couplings. 

We want to stress that our work is yet another indication 
that separating the contributions of the backbone and its envi- 
ronment provides relevant and non-trivial information about 
the nature of the critical behavior of the ±J EA model. We 
believe that this separation should also be important in the 
study of other physical quantities. Moreover, this analysis can 
be applied to the other systems (e.g. K-SAT problem 12411 . 
Viana-Bray spin glass B25I1 . etc.) that are known to have a 
backbone. 

We are working to extend the concept of backbone to sys- 
tems with nondegenerated GSs, as happens with continuous 
distributions of bonds. Results are still too preliminary to be 
reported here. 
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